Location: BondGraph Basic @ 994dbbb88a76 / BG Tutorial Biochemical Systems / TMDD.cellml

Author:
Soroush <ssaf006@aucklanduni.ac.nz>
Date:
2021-04-08 15:18:19+12:00
Desc:
modified svg
Permanent Source URI:
https://models.fieldml.org/workspace/43b/rawfile/994dbbb88a761fe3a62bfaf8352da94103650bb2/BG Tutorial Biochemical Systems/TMDD.cellml

<?xml version='1.0'?>
<model name="TMDD" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#">
    <units name="day">
        <unit multiplier="21600" units="second"/>
    </units>
    <units name="l_per_day">
        <unit units="litre"/>
        <unit exponent="-1" units="day"/>
    </units>
    <units name="nmol_per_l">
        <unit prefix="nano" units="mole"/>
        <unit exponent="-1" units="litre"/>
    </units>
    <units name="per_nmol_per_l_day">
        <unit exponent="-1" units="nmol_per_l"/>
        <unit exponent="-1" units="day"/>
    </units>
    <units name="nmol_per_l_day">
        <unit units="nmol_per_l"/>
        <unit exponent="-1" units="day"/>
    </units>
    <units name="nmol">
        <unit prefix="nano" units="mole"/>
    </units>
    <units name="per_day">
        <unit exponent="-1" units="day"/>
    </units>
    <component name="state">
        <variable name="t" units="day"/>
        <!-- State variables-->
        <variable name="C_tot" units="nmol_per_l"/>
        <variable name="R_tot" units="nmol_per_l"/>
        <variable initial_value="0.0" name="A_d" units="nmol"/>
        <variable initial_value="0.0" name="A_t" units="nmol"/>
        <variable initial_value="33.3" name="C" units="nmol_per_l"/>
        <variable initial_value="5" name="R" units="nmol_per_l"/>
        <variable initial_value="0.0" name="RC" units="nmol_per_l"/>
        <!-- Parameters-->
        <variable initial_value="100" name="Dose" units="nmol"/>
        <variable initial_value="0.45" name="Q" units="l_per_day"/>
        <variable initial_value="0.15" name="CL" units="l_per_day"/>
        <variable initial_value="0.6" name="F_SC" units="dimensionless"/>
        <variable initial_value="1.0" name="k_abs" units="per_day"/>
        <variable initial_value="8.0" name="k_on" units="per_nmol_per_l_day"/>
        <variable initial_value="8.0" name="k_off" units="per_day"/>
        <variable initial_value="1" name="k_syn" units="nmol_per_l_day"/>
        <variable initial_value="0.2" name="k_deg" units="per_day"/>
        <variable initial_value="0.04" name="k_int" units="per_day"/>
        <variable initial_value="3.0" name="V_c" units="litre"/>
        <variable initial_value="3.0" name="V_p" units="litre"/>
        <variable name="k_elim" units="per_day"/>
        <variable name="k_pt" units="per_day"/>
        <variable name="k_tp" units="per_day"/>
        <!-- State equations-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>A_d</ci>
                </apply>
                <apply>
                    <times/>
                    <apply>
                        <minus/>
                        <ci>k_abs</ci>
                    </apply>
                    <ci>A_d</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>A_t</ci>
                </apply>
                <apply>
                    <minus/>
                    <apply>
                        <times/>
                        <ci>k_pt</ci>
                        <ci>C</ci>
                        <ci>V_c</ci>
                    </apply>
                    <apply>
                        <times/>
                        <ci>k_tp</ci>
                        <ci>A_t</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>R</ci>
                </apply>
                <apply>
                    <plus/>
                    <apply>
                        <minus/>
                        <apply>
                            <minus/>
                            <ci>k_syn</ci>
                            <apply>
                                <times/>
                                <ci>k_deg</ci>
                                <ci>R</ci>
                            </apply>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_on</ci>
                            <ci>C</ci>
                            <ci>R</ci>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <ci>k_off</ci>
                        <ci>RC</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>RC</ci>
                </apply>
                <apply>
                    <minus/>
                    <apply>
                        <times/>
                        <ci>k_on</ci>
                        <ci>C</ci>
                        <ci>R</ci>
                    </apply>
                    <apply>
                        <times/>
                        <apply>
                            <plus/>
                            <ci>k_int</ci>
                            <ci>k_off</ci>
                        </apply>
                        <ci>RC</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>C</ci>
                </apply>
                <apply>
                    <plus/>
                    <apply>
                        <minus/>
                        <apply>
                            <minus/>
                            <apply>
                                <divide/>
                                <apply>
                                    <times/>
                                    <ci>k_abs</ci>
                                    <ci>A_d</ci>
                                </apply>
                                <ci>V_c</ci>
                            </apply>
                            <apply>
                                <times/>
                                <apply>
                                    <plus/>
                                    <ci>k_elim</ci>
                                    <ci>k_pt</ci>
                                </apply>
                                <ci>C</ci>
                            </apply>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_on</ci>
                            <ci>C</ci>
                            <ci>R</ci>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <ci>k_off</ci>
                        <ci>RC</ci>
                    </apply>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <ci>k_tp</ci>
                            <ci>A_t</ci>
                        </apply>
                        <ci>V_c</ci>
                    </apply>
                </apply>
            </apply>
            <!-- Parameters equations-->
            <apply>
                <eq/>
                <ci>k_elim</ci>
                <apply>
                    <divide/>
                    <ci>CL</ci>
                    <ci>V_c</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>k_tp</ci>
                <apply>
                    <divide/>
                    <ci>Q</ci>
                    <ci>V_p</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>k_pt</ci>
                <apply>
                    <divide/>
                    <ci>Q</ci>
                    <ci>V_c</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>C_tot</ci>
                <apply>
                    <plus/>
                    <ci>C</ci>
                    <ci>RC</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>R_tot</ci>
                <apply>
                    <plus/>
                    <ci>R</ci>
                    <ci>RC</ci>
                </apply>
            </apply>
        </math>
    </component>
</model>